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We report measurements of the branching fractions and CP asymmetries for B — > wh ± (h = 
K or 7r) and the observation of the decay B° — > rjK° from the final data sample of 772 x 10 6 
BB pairs collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. The 
measured branching fractions are B(B ± nK^ = (2.12 ± 0.23 ± 0.11) x 10~ 6 , B(B ± -> r?vr ± ) = 
(4.07 ± 0.26 ± 0.21) x 10~ 6 and B(B° -» nK°) = (1.27±g;|| ± 0.08) x 10" 6 , where the last decay 
is observed for the first time with a significance of 5.4 standard deviations (<r). We also find 
evidence for CP violation in the charged B modes, Acp(B ± — > nK ± ) = —0.38 ± 0.11 ± 0.01 and 
Acp(B ± — > = —0.19 ± 0.06 ± 0.01 with significances of 3.8u and 3. Oct, respectively. For all 

measurements, the first and second uncertainties are statistical and systematic, respectively. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charge-Parity (CP) violation plays an important role 
in any explanation of the observed dominance of mat- 
ter over antimatter in our Universe. Current experimen- 
tal knowledge about CP violation is limited. Charmless 
hadronic B decays constitute sensitive probes for CP vi- 
olation in the standard model (SM) as well as beyond, 
and can help to further elucidate this unsolved question. 
In the SM, the decays B ± -> jjWjfi and B° -> r)^K° 
are expected to primarily proceed through b — >• s pen- 
guin processes and a b — > u tree transition as shown in 
Fig. Q] The large B -> rfK^% and small B -> r)K%% 
branching fractions can be explained by ?; — r\' mixing 
along with constructive and destructive interference be- 
tween the amplitudes of the two penguin processes 

The branching fraction of B° r]K° is expected to 
be lower than that of B ± — »• r/K^, because the tree di- 
agram in the B° — > r/K° decay is color suppressed. The 
destructive combination of penguin amplitudes may in- 
terfere with the tree amplitude in B — > nK , resulting in 
a large direct CP asymmetry (Acp) [E @, defined as 



A 



CP 



T[B-(B°) -» 7?/t-W] - T[B + (B°) ->■ yyfe+W] 
T[B-(B°) -> rjhFW] +T[B+{B°) -»• rjh+W] ' 



(1) 

where T(B — > rjh) is the partial width obtained for the 
B — > rjh decay, and h denotes K or n. Similarly, direct 
CP violation could be sizeable for B ± — > rjir^ owing to 
the interference between b — > d penguin and b — > u tree 
diagrams. Several theoretical calculations with different 
mechanisms [7- 12 suggest a large Acp for both B — > r/K 
and B — > tjtt, although the sign could be either positive or 
negative. Previous Belle [H and BaBar Q measurements 



indicate a large negative Acp in B^- — > r/K^, but more 
data are needed to be statistically sensitive to a non-zero 
A C p in B ± -> tjtt ± . 

In this Letter, we report the first observation of B° — > 
rjK , and evidence for direct CP asymmetries in B — > 
r/K^ and B ± — ► 777r ± using the final Belle data set. The 
data sample corresponds to (772 ± 11) x 10 6 BB pairs 
collected with the Belle detector at the KEKB e + e~ 
asymmetric-energy (3.5 GcV on 8.0 GeV) collider 



13] 



operating at the T(45) resonance. The production rates 
of B + B~ and B°B° pairs are assumed to be equal at the 
T(45) resonance. 

The Belle detector [14j is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector, a 
50-layer central drift chamber (CDC), an array of aero- 
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FIG. 1: (a) b — > u tree diagram for B + — ¥ r/'' K + . (b) color- 
suppressed b — > u tree diagram for B° — >• n ( ^K°. (c), (d) 
b — > s gluonic penguin diagrams for B -> rj^K. 
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gel threshold Chcrenkov counters (ACC), a barrel-like ar- 
rangement of timc-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprising CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux- return 
located outside the coil is instrumented to detect K\ 
mesons and to identify muons. 

The event selection and B candidate reconstruction 
methods are similar to those described in Ref. We 
select 77 and ir° candidates through the decay chains 
77 — > 77 (f? 77 ), 77 — > Tr + ir~ir° (773^) and 7r° — > 77. 
We require the two photons from the 77 and tt° decays 
to have energies (Eyi, i = 1,2) greater than 50 MeV. 
Candidate ir° mesons are selected from pairs of pho- 
tons with invariant masses between 115 MeV/c 2 and 152 
MeV/c 2 . In the ?y 77 reconstruction, the photon energy 
asymmetry |i? 7 i — i? 7 2|/(.E 7 i +£^,2), is required to be less 
than 0.9 to reduce the large combinatorial background 
from low-energy photons. Photons in 7y 77 reconstruction 
are not allowed to pair with any other photon having 
Ej > 100 MeV, to form a n° candidate. We require the 
invariant mass of the 77 77 and 773^ candidates to be in the 
intervals (501,573) MeV/c 2 and (538.5,556.5) MeV/c 2 , 
respectively. In order to improve the 7r° and 77 energy 
resolution, a mass-constrained kinematic fit is performed 
after the candidate selection. 

Charged tracks that are not used to form K s can- 
didates (see below) are required to have a distance of 
closest approach with respect to the interaction point 
(IP) of less than 3.0 cm along the electron beam direc- 
tion (z) and less than 0.3 cm in the transverse plane. 
Charged kaons and pions are identified using dE/dx 
information from the CDC, Cherenkov light yields in 
the ACC, and time-of-flight information from the TOF. 
This information is combined to form a likelihood ratio, 
T^k/tt = C-k/{£k + £tt), where Ck (At) is the likeli- 
hood of the track being a kaon (pion). Charged tracks 
with TZk/w > 0.6 (< 0.4) are treated as kaons (pions) for 
£? ± — > r\K ± [B ± — > ?/7r ± ) reconstruction. A less strin- 
gent requirement, IZk/tt < 0-6, is used for charged pions 
in the 773^ selection. The efficiencies for the TZk/tv rG ~ 
quirement are 84% for kaons, 89% for pions, and 94% for 
pions in 773^ reconstruction. Furthermore, for B — > rjh^ 
and B° — >• rjK° we reject charged tracks consistent with 
either the electron or muon hypothesis. 

Candidate K° mesons are reconstructed in K% —> 
7r + 7r~ decays. The Kg candidates are required to have 
an invariant mass lying between 488 MeV/c 2 and 508 
MeV/c 2 . The charged tracks for each Kg candidate are 
required to have a distance-of-closest approach with re- 
spect to the IP of larger than 0.02 cm in the transverse 
plane. The angle between the K s momentum and the 
direction from the IP to the K s decay vertex must be 
within 0.03 rad. The distance between the two daughter 
tracks at their point of closest approach in the transverse 
plane is required to be less than 2.40 cm, and the flight 



length of the Kg is required to be larger than 0.22 cm. 

Candidate B mesons are identified using the 
modified beam-energy-constrained mass -Mb c = 
\/(^bea m /c 2 ) 2 " \P* B /c\ 2 IH, and the energy differ- 
ence AE = E* - £ b * cam , where £ b * cam is the beam 
energy, and Eg and p B are the energy and modified 
momentum, respectively, of the B candidate in the 
T(45) rest frame. The energy Eg is calculated as 
E* B = E* + E%. The momentum p B is calculated 
according to 

P*b = P* h + x y/(E£ cam E* f - ml (2) 

where m, is the nominal 77 mass [lij . Since charged 
tracks are generally measured with a better precision 
than photons, the 77 decays to neutral particles have worse 
momentum resolution than primary charged tracks from 
B decays. The p* B resolution is improved using Eq. (2) 
because p" h and E£ eam are determined more precisely than 
p* v . Events with M bc > 5.2 GcV/c 2 and \AE\ < 0.3 GeV 
are retained for further analysis. 

The dominant background arises from e + e~ — > qq (q = 
u, d, s, c) continuum events. We use event topology vari- 
ables to distinguish spherically distributed BB events 
from the jet-like continuum background. First we com- 
bine a set of modified Fox- Wolfram moments [ItJ into a 
Fisher discriminant. We then compute a likelihood from 
the product of probability density functions (PDFs) that 
describe the Fisher discriminant, cos#J, and Az distri- 
butions. Here, 0* B is the angle between the B flight di- 
rection and the beam direction in the T(4S) rest frame, 
and Az is the decay flight-length difference, along the z 
axis, between vertices of the signal B and the accompany- 
ing B. A likelihood ratio, 1Z = C S /(C S + C q q), is formed 
from signal (C s ) and background { Cqg ) likelihoods, which 
are obtained from GEANT-based Monte Carlo (MC) 
simulated samples. Signal MC events are generated with 
EVTGEN 0, which invokes the PHOTOS package 
to take final state radiation into account. We require 
1Z > 0.2 to suppress continuum background in all modes 
and then translate 71 to 1Z' , defined as 

U' = In { ^~ n ^l ) ■ (3) 

In this expression TZ m \ n (TZ ma x) is equal to 0.2 (1.0). This 
translation is convenient, as the TZ' distributions for sig- 
nal and backgrounds can be described by a simple sum 
of Gaussian functions. 

Signal yields are extracted by performing an unbinned 
extended three-dimensional maximum likelihood fit. The 
likelihood for each B + mode is defined as 

C = < v ■ II X V 7 ' 

i 3 

V) = \[l-q l A CP] ]P 3 {AP hcl AE\1Z' 1 ), (4) 
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TABLE I: Detection efficiency (e) including sub-decay branching fractions, yield, significance of the yield 2(3*0, measured 
branching fraction B, charge asymmetry Acp, and significance of the charge asymmetry Yi(Acp) for B — > r\h decays. The first 
errors are statistical and the second ones are systematic. 
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FIG. 2: A 75 (left), M bc (middle) and ft' (right) distributions for 73° ->■ t^Ts" candidate events with the 77 77 and 773^ modes 
combined. Points with error bars represent the data, the total fit functions are shown by black solid curves, signals are shown 
by red solid curves, dashed curves show the continuum contributions and filled histograms are the contributions from charmless 
B backgrounds. The AE, Mb c , and ft' projections of the fit are for events that have 5.27 GeV/c 2 < Mbc < 5.30 GeV/c 2 and 
ft' >0.55, -0.10 GeV< AE < 0.08 GeV and ft' >0.55, and -0.10 GeV< AE < 0.08 GeV and 5.27 GeV/c 2 < M bc < 5.30 
GeV/c 2 , respectively. 



where i denotes the i event. Category j runs over the 
signal and background components, where background 
components include continuum, the cross-feed due to 
K-ir rnisidcntification, and the background from other 
charmless B decays, but do not include b — > c decays, 
which arc found to give a negligible contribution. The pa- 
rameter Nj represents the number of events for category 
j, Pj(Mi c> AE i ,TZ' i ) is the PDF in M hc , AE and TZ', and 
q is the i?-meson charge. For the 73° mode, the Acpj pa- 
rameters are set to zero since CP asymmetries cannot be 
determined without additional information. The validity 
of the three-dimensional fit is tested with large ensemble 
tests using MC samples and by fits to a high-statistics 
control sample of B + D°n + (D° -> K + n~iT ) decays. 

All the signal and cross-feed background PDFs in M^ c 
and TZ' are modeled with a single Gaussian. In 73 — > r)jyh 
(B —> r)3 T h) modes, the PDFs in AE are described by a 
Crystal Ball [2l[ (a sum of two Gaussians) function . The 
peak positions and resolutions in M^ c , AE and TZ' are 
adjusted according to the data-MC differences observed 
in control samples [HI- 

The continuum background in AE is described by 
a second-order polynomial, while the M\ JC distribution 
is parameterized with an ARGUS function, f(x) = 



xy/1 — x 2 exp [— £(1 — x 2 )}, where x is Mb c /Eb eam and £ 
is a free parameter in the fit (23j. The TV PDF is a dou- 
ble Gaussian function. The background PDFs in both 
Mbc and AE for charmless B decays are modeled with 
smoothed two-dimensional histograms obtained from a 
large MC sample, and the TV PDF is a single Gaussian. 

We perform a simultaneous fit to 73+- — > r\K^ and 
73+- — > rjir^ candidate events, since these two decay 
modes can feed into each other. In the likelihood fits, 
Nj and Ac pj arc allowed to vary for the continuum and 
charmless B backgrounds. Shape parameters of the con- 
tinuum background PDFs are also floated. The cross-feed 
background in r\K ± (rjn^) and signal in r/n^ (rjK^) share 
the same fitting parameters in both Aqp and branching 
fraction. Figure [2] shows the Mb c , AE 1 and TZ' projections 
of the fit to the B° -> r]K° s sample. The M hc and AE 
projections for the 73+ — > ??/i+ and 73~ — > r]h~ samples 
are shown separately in Fig. [3] 

The branching fraction for each mode is calculated 
by dividing the efficiency-corrected signal yield by the 
number of -B73 pairs. The dominant systematic er- 
rors on the branching fraction come from MC model- 
ing of the rj, 7T°, and Kg selection efficiency; these er- 
rors are 4.0%, 4.0%, and 1.6%, respectively. The sys- 
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tematic error due to TZ(K/%) requirement is 0.9% for 
kaons and 0.8% for pions. It is estimated from the 



") sample. The systematic 



D*+ -> D°n+(D° 
error due to the charged-track reconstruction efficiency 
is estimated to be 0.35% per track, which is determined 
from a study of the D*± -> D° 7 r ± (Li -> tt+tt"^) de- 
cay. Any difference in the efficiency when the 1Z crite- 
rion is applied to data or MC is investigated using the 
B+ -> D°ir+(D° -> A'+tt-tt ) sample; the results of 
this study imply a 0.6% systematic uncertainty. The fit- 
ting systematic errors due to the signal PDF modeling 
are estimated from changes of the fit parameters while 
varying the calibration factors by one standard devia- 
tion. The systematic error due to charmless B back- 
ground PDF modeling is calculated from the difference 
observed between the signal yield when the charmless 
yield is floated in the fit and that when the yield is fixed 
to the MC expectation. The systematic error due to the 
uncertainty in the total number of BB pairs is 1.4%, 
and the error due to limited signal MC statistics used 
to evaluate the efficiency is 0.55%. The systematic er- 
rors on Aqp arise from detector bias, uncertainties on 
the detector bias and PDF modeling. The possible de- 
tector bias due to the tracking acceptance and lZ(K/ir) 
selection for Acp{B ± —> tjtt ) is evaluated u sing the 
fitted Aqp value of the continuum background |24l |25| . 
The detector bias in Aqp{B ± — > -qK^) is evaluated us- 
ing the D+ -> 4)ir+(4> -> K+K~) and D° -> K~-k+ 



samples [2J, |25j. There is a contribution to the Acp 
systematic uncertainty from the modeling of the signal 
PDFs. The total systematic errors for Aqp are in the 
range (8.2 - 14.2) x 10" 3 . 

The statistical significance is evaluated as 
yj — 2 In ( Co / £ max ) i where Co is the likelihood value 
when either the signal yield or Aqp is fixed to zero, 



and C T 



is the nominal likelihood value. The total 



significance (S) including PDF modeling systematic 
uncertainty is calculated after smearing the likelihood 
distribution with the appropriate PDF modeling sys- 
tematic error. In Table Q] we list the fitted signal yields, 
charge asymmetries, reconstruction efficiencies, and 
branching fractions. The combined result for the two rj 
decay modes is obtained from the combined likelihood 
function. 

In summary, using the final Belle data sample con- 
taining 772 x 10 6 BB pairs and a three-dimensional 
fit that maximizes the efficiency, we provide new mea- 
surements based on signal yields 2.5 times larger than 
those reported in our previous publications (jj . We 
find evidence for CP asymmetries in B^ — > rjK^ and 
B ± -> t]^: A CP (B ± -> rjK*) = -0.38 ± 0.11 ± 0.01 



and A CP (B ± -> rpr*) = -0.19 ± 0.06 ± 0.01. The sig- 



nificance of Acp(r]K + ) [Acp(r]ir + )] is 3.8a [3. Oct] 



Ev- 



idence for Acp(B ± — > fj'K ) is seen for the first time. 
We also observe the decay B° 7]K° for the first 
time with a significance of 5.4ct and a branching fraction 



-B + -> V K+ | 
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FIG. 3: AA and M hc projections for B + -> r]h + (left) and 
B~ — > r;h~ (right) candidate events with the ^ 77 and 773^ 
modes combined. Points with error bars represent the data, 
the total fit functions are shown by black solid curves, sig- 
nals are shown by red solid curves, dashed curves are the 
continuum contributions, dotted curves are the cross-feed 
backgrounds from misidentification and filled histograms are 
the contributions from charmless B backgrounds. The AA 
and Mbc projections of the fits are for events that have 
5.27 GeV/c 2 < M hc < 5.30 GeV/c 2 and K' >1.95, -0.10 
GeV< AA < 0.08 GeV and 11' >1.95, respectively. 



B(B° -+ r,K°) = (1.2711$; ± 



x 10" 



In addition, 



we report the following new measurements of the branch- 
ing fractions: B(B ± -> t/K^) = (2.12±0.23±0.11)xl0- 6 
and B(B ± -^■riir ± ) = (4.07±0.26±0.21) x 10" 6 . All our 
branching fraction and Aqp measurements supersede the 
results in Ref.0]. 
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